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Abstract. GRBs are nowadays a rather well understood phe- 
nomenon in the soft (KeV-MeV) 7-ray energy band, while 
only a few GRBs have been observed at high photon energies 
{E~f > 1 GeV). It is also widely recognized that GRBs acceler- 
ate protons to relativistic energies and that dense media are 
often present nearby the sources. Within this framework and 
by further adopting Totani's suggestion that GRB events re- 
lease an amount of energy ~ 10^*Af2 erg, we compute in detail 
the high-energy 7-ray flux from the decay of neutral pions pro- 
duced through the interaction of accelerate protons with nu- 
cleons in the surrounding medium. We also take into account 
the local and intergalactic 7-ray absorption. The presence of 
magnetic fields around the GRB sources causes the deflection 
of the accelerated protons and so a temporal spread of the 
produced high-energy 7-rays with respect to the signal in the 
soft 7-ray band. Moreover, we analyze the possibility to detect 
the 7-ray signal in the GeV-TeV energy range by the ARGO 
detector under construction in Tibet. 
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1. Introduction 

Gamma-ray bursts (GRBs) have intrigued observers and con- 
founded theorists ever since their discovery over thirty years 
ago. Up to now, the most meaningful contribution to the com- 
prehension of the GRB physics has been obtained through the 
observations by BATSE and EGRET instruments on board 
of the CGRO satellite. Since 1991, BATSE detects about 
one GRB per day in the 50 - 300 keV energy band (here- 
after indicated as soft 7-ray regime), with fluences between 
~ 10~^ — 10~^ erg cm~^ a nd du rations ranging from ~ 10 ms 
to ~ 10^ s (Paciesas et al. |l999| ). 

The accumulation of a substantial GRB population by 
BATSE shows that they are isotropically distributed. The sim- 
plest explanation for this fact is the GRB cosmological origin, 
at variance with the the pre-CGRO view in which GRBs were 
associated to neutron stars distributed in a thick galactic disk. 

The GRB cosmological origin has been clamorously con- 
firme d mor e recently by the Beppo-SAX observations (Costa 
et al. 1997) of counterparts in the X-ray band of some GRBs. 
Thanks to these observations, it has been possible not only a 
X-ray identification of GRBs, but also the determination with 



a better accuracy of the GRB source position, allowing the 
conse quent tracking of optical and radio telescopes (Galama 



et al. 



1997) 



Table 1. For some bursts for which it has been possible to 
estimate the distance, we give their redshift z and the energy 
5^"° radiated in the soft 7-ray band. In the third column we 
indicate how the GRB redshift has been measured: (a) from the 
redshift of the host galaxy; (b) from Lyman break of the optical 
afterglow; (c) from absorption lines in the optical afterglow. 



GRB 



970228 
970508 
971214 
980329 
980613 
980703 
990123 

990510 
991208 
991216 



0.695 
0.835 

3.42 

~5 
1.0964 
0.966 
>1.6 

1.62 
0.7055 

1.02 



Note 



5.2 X lO"" 
~ 7 X 10^^ 
~ 3 X 10^^ 
>10^'' 

5.2 X 10=^ 
>10^^ 

> 3 X 10^^ 

2.9 X 10^^ 

1.03 X 10^^ 
6.7 X lO'^^ 



Ref. 



Djorgovski |l999a 
Metzger 
Kulkarni 
Fruchter 

Djorgovski [19991: 
Djorgovski 199< 
Kulkarni |l999 




Harrison 199S| 



Sagar 
Sagar 



200c 



200c 



In turn, redshift measurements (through optical observa- 
tions) allow to evaluate the distance of some GRBs and there- 
fore to derive the amount of energy they should have radiated 
in an isotropic explosion (henceforth indicated as f^"") in the 
soft 7-ray energy band. 

Currently, the most popular explanation for the G RB phe- 
nomenon is the fireball model (Rees & Meszaros 1992), i. e. the 
dissipation of the kinetic energy of a relativistic (with Lorentz 
factor r ~ 10^ — 10^) expanding blast wave in internal shocks 



(see Piran 1999 for a recent review). Due to the short time 
scale variability of GRBs (> 10""^ s), soft 7-rays are generally 
considered to be emitted by relativistic electrons, through syn- 
chrotron and inverse Compton mechanisms. This model is par- 
ticularly successful in explaining the long wavelength afterglow 
behaviour observed for some GRBs. 

In Table 1, for some bursts for which it has been possible to 
estimate the distance, we give their redshift z and the inferred 
f "° values. Note that if GRBs are beamed s ource s as indicated 
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by afterglow observations (Huang, Dai & Lu 199E| ) , the emitted 
energy would be decreased by the factor AQ/Att, where Af2 is 
the beaming angle. As one can see in Table 1, the brightest 
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burst (GRB 990123) has f^"" > 3 x 10^'' erg showing that in 
a GRB event up to one solar mass can be radiated in form of 
soft 7-rays. 

It is also well known that GRBs are sources of high-energy 
7-rays, as it has been shown by the EGRET spark chamber, 
which detected 7-rays in the energy band ~ 30 MeV - 20 GeV, 
from a high percentage of bri ght B ATSE bursts (Schneid et al. 



1992 



1995 and Hurley et al. 1994). These observations clearly 



show that the GRB phenomenon is not exclusively the domain 
of the soft 7-ray energy band. 

The natural question which raises is the production mech- 
anism of the high-energy 7-rays (with energy > 1 GeV) 
that, in principle, may be due to both electrons and protons. 

However, if the electron emission spectra which fit BATSE 
observations are extrapolated in the high-energy domain, the 
obtained photon fluxes at Earth a re we ll below the values de- 
tected by EGRET (PiUa & Loeb |199^). This fact manifestly 



implies that high-energy 7-rays from GRBs cannot be entirely 
accounted for by radiation emitted by accelerated electrons. 

Indeed, as is usual for many high-energy astrophysical 
sources (such as AGN, Seyfert galaxies, etc.), also for GRBs 
it is expected that the most efficient mechanism to produce 
high-energy 7-rays is via n" decay of pious primarily produced 
through proton-nucleon (pN) interactions. EJ 

Within the framework of the fireball model for GRBs, it 
has been recently argued that all GRB events could release 
roushl¥| the same amount of energy £ ~ 10^'^ AQ erg (Totani 



199£) 



The total GRB energy is shared between electrons and pro- 
tons which, at least in the initial stage of internal shock gen- 
eration, carry a much larger amount of energy than electrons, 
by a factor mp/m^ ~ 2000. 

It is uncertain what fraction of the initial energy is subse- 
quently transferred by protons to electrons, but the simplest 
Coulomb interaction cannot be very efficient within the short- 
est GRB time-scale variability of ~ 10~^ s. Another mecha- 
nism to transfer energy to electrons could be through proton 
synchrotron emission and subsequent annihilation 77 — > e^e~, 
but also this process is generally inefficient owing to-iits quite 
long time-scale compared with the GRB variability, u 



^In principle, accelerated protons could radiate their energy 
through synchrotron emission, before strong interactions are in 
operation. However, synchrotron emission is relevant o nlv for 



Indeed, according to the fireball model, the energy transfer 
from the proton into the electron component occurs via internal 
shock generation. The efficiency of this mechanism depends 
sensitively on the fireball Lorentz factor F so that the spread 
of the r values should account for the spread in the observed 
soft 7-ray emissivity E-^ of the electronic component. 

However, as one can see from Table 1 at least for the few 
GRBs for which it has been possible to estimate the distance, 
the energy is a minor fraction (between ~ 10 ^ and ~ 0.1) 



of the total GRB energetic assumed by Totani (1999). So, it 
is likely that a substantial fraction of £ remains in form of 
accelerated protons coming out from the GRB source with a 
total energy £p ~ lO^^'Ail erg. This is the reference value that 
we will adopt in our calculations. 

It is by now widely recognised that GRBs ac celera te pro- 
tons at energies even up ^ 1 0^° eV (Waxman 1995| , Vietri 
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Bottcher & Dermer 1998). If it is so, these protons could 
in turn interact with nucleons in the medium surrounding the 
GRB source, giving rise to tt" and ultimately to a high-energy 
7-ray fiux on Earth. Clearly, a key requirement of this model 
is the existence of a dense enough cloud (which serves as a 
target for pN interactions) surrounding the GRB source or, 
alternatively, the existence of a cloud along the line of sight 
observer-burst and close enough to the source. Moreover, as 
we shall see in the following Sections, in order to have an ob- 
servable high-energy 7-ray flux on Earth, the cloud number 
density njv has to be in the range 10* - 10" cm . 

At flrst sight, so high densities may seem unlikely. How- 
ever, the presence of high density regions nearby GR B sou rces 



has been invoked in several cases. For example, Katz 1994 has 
shown that a cloud with density njv ~ 2 x 10^^ cm~'^ and 
size ro ~ 10^^ cm is required to explain the delayed high- 
energy photons observed more than an hour following the GRB 
940217. Moreover, Piro et al. ( |l999[ ) and Yoshida et al. ( |l999| ) 
have reported an iron emission line in the X-ray afterglow spec- 
trum of GRB 970508 and of GRB 970828. The obtained line 
intensity requires a medium of density un — 10^" cm~'^, thick- 
ness AR ^ lO^ '* cm and size 10^^ cm (Lazzati, Campana & 
Ghisellini 199£). The involved geometry of the target region 



protons with extremely high-energy ~ 10^° eV (Vietri 1997) ers (Paczynski 



could be either a cloud with a significant covering factor lo- 
cated off the line of sight burst-observer or a homogeneous 
spherical shell centered in the GRB progenitor. 

Three main classes of models have been proposed for the 
origin of GR Bs: ne utron star-neutro ri sta rs (NS-NS) merg- 



1986 



and Eichler et al. 



and it is negligible for the production of ~ 1 GeV - 10 TeV 
photons, which are produced much more efficiently through 
pN interactions. 

^We would like to mention that Kulkarni et al. (2000) have sug- 
gested that afterglow observations may offer a robust method 
to evaluate the GRB energetics. In particular, from late time 
X-ray observations and assuming that the fraction of shock en- 
ergy carrie d by electrons is close to equipartition, Freedman & 
Waxman ( |l99g| ) estimate the GRB energy release £ to to be 



1989), "hypernov ae" or 



failed type lb supernovae (Woosley 1993| and Paczynski 1998) 



in the range 3 x 10 — 3 x 10 erg, which must be considered 
a more conservati ve ra nge of (low) GRB energetics. Ho wever , 
both the Totani (|l999|) and the Freedman & Waxman ([1999D 



and "supranovae" (Vietri & Stella |l998| ). In the NS-NS model 
the target cloud could be associated to an accretion disk around 
the GRB source. In the hypernova case, the burst is caused by 
the evolution of a massive (~ 100 Mq) star located in a dense 
molecular cloud, while in the supranova scenario, a shell-like 
supernova remnant is naturally left over around the location 
of the GRB. 

In all the above mentioned scenarios, high-energy 7-ray 
production via pN interactions is expected with an inten- 
sity roughly proportional to un. Moreover, a prediction of 
our model is that a long wavelength afterglow emission from 



estimates for £ may still be considered consistent with obser- 
vations, due to the uncertainties and precarious nature of GRB 
energy estimates (see also Kulkarni et al. ^000| ) . 
■^It has been shown that the proton cooling time becomes com- 
parable to the GRB duration only at the maximum proton 



energy (~ 10"^" eV), and, even in this case, the energy radi- 
ated in the GeV- TeV range is considerable only if the proton 
energy spectrum is harder than th e tvpical shoc k acceleration 
spectrum Jp{Ep) oc E'^-'^ (Totani |l9984 |l998b| ). 
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GRBs should be present only if the cloud covering factor is low 
enough. 

A feature of the model in question is a delay in the arrival 
times of the high-energy 7-rays with respect to the soft 7- 
rays and a temporal spread of the signal, as a consequence of 
the proton deflection due to the presence of magnetic fields in 
the region surrounding the GRB source. A lower limit to the 
magnetic field is surely the typical interstellar magnetic field 
BisM — 1 /J-G. 

The aim of the paper is to estimate in detail the flux on 
Earth of 7-rays in the GeV-TeV energy range from cosmologi- 
cal GRBs and in particular the temporal structure of the signal. 
In Sections 2 and 3 we present the model we use to calculate 
the source function for 7-ray production via pN interactions. 
In Section 4 we present our results for the 7-ray flux on Earth, 
taking into account the intergalactic 7-ray absorption. In Sec- 
tion 5 we briefly discuss the ground-based high-energy 7-ray 
experiments in connection with GRB observations and give the 
model parameter range for which the GRB high-energy detec- 
tio n may be possible by the ARGO-YBJ experiment (Abbres- 
Finally, our conclusions are presented in 



199^ 



cia 1996, Bacci 
Section 6. 

2. The Model 



As previously stated, we assume that the GRB source emits an 
amount of energy erg in the form of relativistic 

protons, released during a time At ~ li-e, at a distance ro ~ 
10^^ cm from the central GRB source We further assume 
that there exists a dense enough cloud near/around the GRB 
source so that accelerated protons subsequently interact with 
nucleons in the cloud giving rise to tt" and ultimately to high- 
energy 7-rays. 

An essential ingredient to calculate the 7-ray production 
rate is the proton flux Jp{Ep,r) which is a function of both 
the proton energy and the radial coordinate r. As impli ed by 



Fermi type shoc k acc eleration models (see e.g. Gaisser I990|, 



Berezinskii et al. 1990), Jp{Ep, r) is assumed to follow a power- 
law energy spectrum. It moreover depends on the distance r 
from the central engine according to a function /(r) to be later 
specified 



Jp{Ep,r) = CiEp/Gey)--''f(r) 

protons cm~^ s~^ sr~^ GeV^'^ 



(1) 



with ultrarelativistic velocities and further diffusive processes 
are negligible during the time before pN interactions, w e will 
adopt the r elayant value Up ~ 2.2 (Bednarz & Ostrowski 1998, 
Vietri |2000| ).a 



The radial profile /(r) depends of course on the geometrical 
dilution factor which follows the dependence. Moreover, we 
take into account that protons emerging from the GRB source 
at the distance ro, subsequently contribute to the production 
of high-energy 7-rays just in their first pN interaction (for sim- 
plicity we neglect interactions of secondary particles with lower 
energy). The probability that this first interaction occurs af- 
ter a path of length I is equal to exp(— cTpjvnjvZ), where apjv 
is the total pN cross-section (o-pjv ~ 3 x 10~^^ cm'^, nearly 
independent of the proton energy) and njv is the nucleon num- 
ber density of the cloud. Correspondingly, the average length 
travelled by protons in the surrounding medium results to be 
AR = (o-pivnjv)"^ ~ 10"^ pc (10^° cm"^/niv). The mass en- 
closed wiU be A/ci ~ 0.4 MQ{ro/10^^ cmf in a shell-like cloud 
geometry and Md ~ 0.2 Mq{W^'^ cm^'^/wiv)^ in the case of 
an intervening cloud. 

An important point to stress is that high-energy 7-rays 
may in principle be re-absorbed if the path they have to travel 
inside the cloud (after being produced) is long enough. If the 
dominant absorption process is through -yH interactions, the 
corres pondin g column density is ~ 3x 10^^ cm~^ (Particle Data 
Group 1996 ) or, equivalently, an attenuation length I^h ~ AR. 
In this case high-energy 7-rays should suffer attenuation. How- 
ever, it is expected that clouds nearby the GRB source present 
a high ionization degree, so that the relevant absorption pro- 
cesses are 76 and yp interactions, which have a much smaller 
cross-section implying a large 7 free-path before re-absorption. 

Another effect we have to consider in defining the func- 
tion /(r) is the defiection of the proton trajectory before the 
interaction with target nucleons occurs. In fact, protons are 
defiected by the intervening magnetic field B by an angle 
a(Ep) — 1/Rl{Ep), where Rl{Ep) = Ep/{eB) is the Larmor 
radius. From simple geometrical considerations, protons start- 
ing at r = ro and reaching the radial coordinate r are deflected 
by an angle 



a ~ arcsm 



r — ro 



Assuming /(ro) — 1, the constant C is determined by the 
normalization condition which yields 



^Rl{Ep) 
so that / is given by 

l{Ep,r) Rl{Ep) arcsin 



/ r - ro 
\Rl{Ep) 



(4) 



(5) 



C 



Anrl I{ap) At 
where 



I{ap)= / (Ep/GeVy'-EpdEp 

J Eq 



(2) 



(3) 



and -Eo ~ VvripC? is the minimum proton energy. 
As far as the value of the energy spectral index ap is con- 
cerned, since protons are accelerated by shock waves moving 



We have verified that our results do not sensitively depend on 
the adopted values of At and ro. 



Moreover, protons defiected by an angle a > 7r/2, i.e. pro- 
tons that have travelled for a length / > 7r7?_t(£'p)/2, turn back 
towards the GRB source, so they presumably do not contribute 
to the production of observable high-energy 7-rays. 

On the basis of these assumptions, the radial profile takes 
the form: 



In this respect we mention that the observed local CR spec- 
trum has a power index ap ~ 2.75, softer than the previous 
one. As it is well known, this is due to the fact that CRs diffuse 
throughout the Gala:xy, scattering on the hynomogeneities of 
the magnetic field. In this diffusion process, high-energy parti- 
cles rapidly loose energy and at the end of the diffusion process, 
the CR power index gets increased with respect to initial value. 
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/(r) = (J^y exp[-apNnNl{Ep,r)] if / < Rl{Ep) 
/(r)=0 if 1>Rl{Ep) . 



(6) 



The effect of the presence of a magnetic field around the 
GRB source is particularly relevant as regards the 7-ray pro- 
duction rate. The ratio 77 = RL{Ep)crpNnN between the scale 
length for proton deflection by magnetic fields Rl{Ep) and 
that for pN interactions (cTpjvnjv)"^, sets the relative strength 
of the two effects, the deflection and the interaction. 

In the case ij <1, protons can be significantly deflected be- 
fore interacting, so that most of the generated 7-rays will not 
reach Earth. Now, since Rl oc Ep, once um and B are given, 
the condition 77 < 1 means 



Ep < GeV 



10^ 



njv 



B 
1 /iG 



pN -^71^+ anything , tt" 



77 



(8) 



q^{E^f,r) = 16tt UN / dE-^ / dEp Jp{Ep,r)x 

I ^mi-n J 



f 

fJ CO 



dcosO E 



dp*'-^ 



-3 -1 ri 17-1 

cm s GeV 



dp*'' 



= P*\/Pma 

= 6.1, C2 



3.3, 



as reference values, 
)_)_ protons cm~^ s~^ 
ro ~ lO^'* cm, 




(7) 



where the magnetic held at r > ro has been taken of the 
order of magnitude of the interstellar magnetic field. 

3. Source function 

As we discussed in the Introduction the relevant mechanism 
for the high-energy 7-ray production with E~, > 1 GeV is the 
process 



> 

I) 





Following Dermer 198f: , the production spectrum of high- ^ 
energy 7-rays resulting from pN collisions is given by ^ 







-10 



(9) 
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Log(E./'GeV) 
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where 9 is the angle between the line of sight and the 
proton direction and asterisks refer to CM frame. Provided 

-1 < C0s6max < 1, one gets cos 6max = VE-^ ~ E^ax/iPTpTv)- 

Accor ding to Stephens-Badwar's model (Stephens & Bad- 



1981), the Lorentz invariant cross-section for ■k'^ produc- 
tion in pN collisions inferred from high-energy experimental 
data can be parameterized in the form 



(10) 



where x = (xj]^ + (4/s)(p5 
2mp{Ep + mp) and A = 140, B = 5.43, Ci 
C3 = 0.6. 

Numerical values of the 7-ray source function in eq. 
depend on the proton flux Jp{Ep,r) at distance r, the Lorentz 
factor r and the nucleon number density njv. The proton flux, 
in turn, depends on the constant C and the spectral index ap. 

In order to clarify the dependence of q^{E^,r) on the rel- 
evant parameters ap and F we assume 
UN ^ 10^ cm-3 and C ~ W^^/iAnlia 
sr^^ GeV^^ or, equivalently form eq. 
£p ~ lO^^Afi erg and At ~ 1 s. 



Log(E,/GeV) 



Fig. 1. Production spectra qj{E-y) of high-energy 7-rays are given as 
a function of E~f for several values of F and ap, assuming riff = 10^ 
cm-3, and C = lO'^^ /{4:nl{ap)) protons cm'^ s"! sr"! GeV'^ In 
Fig. la, we set F = 1 and give results for ap = 2.2 (full line) and 
2.75 (dotted line). In Fig. lb, for the relevant value ap = 2.2, we see 
the effect on q-y{E'j) of increasing the values of F = 1, 100, 1000. 



In Figs, la and lb we give q^{E^) as a function of E^. In 
Fig. la F is flxed to be one, while ap is 2.2 (continuous line) and 
2.75 (dotted line). In Fig. lb we set ap = 2.2, while F values 
are 1, 100 and 1000. As we can see in Fig. la, consideration of 
a harder proton spectrum increases the high-energy 7-ray pro- 
duction rate (by a factor of ~ 100 at E-y ~ 1 TeV). Moreover, 
we find that the 7-ray spectrum roughly follows a power law 
behaviour qy{Ej,r) ~ S^"^, with power index a^ ~ ap. In 
Fig. lb, we see that an enhanced high-energy 7-ray production 
is also obtained by increasing the F values. 

Production spectra for different values of C and njv can be 
obtained from the results in Fig. 1 by multiplying the quoted 
values by the factor C (47rT(ap))/10^^) x fiiv/lO^. In particular, 
assuming un = 1 cm~"^ and C — 2.1 protons cm~^ s~^ sr~^ 
GeV~^ we can verify that our calculations reproduce, within 
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a few percent, the results on the 7-ray spectrum produced by 
the lo cal co smic-ray interactions with the interstellar medium 
(Mori |1997|). 



4. Gamma-ray flux on Earth 

In order to calculate the differential 7-ray luminosity 
{E^ \ro,R), we have to integrate the source function 
q^{E-^,r) in eq. ^ over the whole volume of the region in 
which the pN interactions occur, i.e. the shell between the 
radii ro and R ~ ro + l/(crpjvnjv): 



Qj{E^;ro,R) ^ q-,{E-,,rydr 7 s^' GeV"^ 



(11) 



Note that here, as stated at the end of Section 2, we have 
neglected re-absorption of the 7-rays in the cloud itself. 

The explicit dependence of Q^{E^\ro, R) on the GRB pa- 
rameters is obtained by using eq. (|l|) in eq. (^), namely 



Q-y{Ej-ro, R) 



47rnjvfp 
I{ap)At 



dEr 



Er, 



f 

•J m 



(12) 



dr ex.p[—apNnNl{Ep,r)] 



dcosO E 



dp*'^ 



As we pointed out in Section 2, the presence of magnetic 
fields around the GRB source affects the differential 7-ray lu- 
minosity in eq. ( |l^ ) through the term l{Ep, r). In Fig. 2, adopt- 




Log(EyGeV) 



Fig. 2. The differential high-energy 7-ray luminosity emitted by a 
GRB source via pN interactions is given as a function of E^, for 
-B = (dotted line) and _B = 1 G (continuous lines). In the lat- 
ter case, the cloud density values are njy = 10*, 10^, 10^°, 10" 
cm~^, from the bottom to the top. Values of Q-y for different 
model parameters, can be obtained by using the scaling relation 
Q^(£;''°/1055 erg)(At/ls)-i. 

ing the same values for the GRB parameters as in Fig. 1 and 
further assuming ap = 2.2, F = 100 and njv in the range 
10** - 10" cm ^, we compare Qj{Ery) obtained in the case 
-B = (in which protons follow straight trajectories) with the 



case 5 = 1 /iG (in which protons of lower energies are de- 
flected). As one can see, in the case B = 0, one obtains the 
same curve for any value of the cloud number density. This 
result is expected since, increasing the density, the scale length 
for the pN interactions decreases, but the projected column 
density Rhn remains constant. 

From Fig. 2 one can see that the effect of the presence of the 
magnetic fields is to diminish the 7-ray production rate, in the 
sense that for a given value of B, Q^(E^\ ro, R) decreases as um 
decreases. In fact, low energy photons are mainly generated by 
pN interactions of low energy protons (generally with energies 
of one order of magnitude above the photon energy) , which are 
more easily deflected by the intervening magnetic flelds and 
therefore give rise to photons that cannot reach the observer. 
There exists, however, a limiting density (~ 10^^ cm~^ in the 
case B ~ 1 ^G) above which most of the protons interact 
before suffer substantial deflection so that the differential 7- 
ray luminosity for B = is recovered. 

Another important effect of the presence of magnetic fields 
is a temporal spread of the emitted high-energy 7-ray signal on 
Earth. To better analyze this point, we can assume a proton in- 
jection time /S.t ~ 1 s, which may be considered instantaneous 
since it is much smaller than the average proton interaction 
time ((jpjvnjvc)"^ ~ 10^ s (10^'' cm~^/njv). Accordingly, de- 
noting by to the instant at which photons arrive on Earth in 
the case B = (no proton deflection), the delay t — to is a 
function of the length I travelled by the protons (before inter- 
acting) and of the deflection angle a given by eq. So we 
get 



r i - ro 



Rl{Ep) / . 

ti - to = arcsm „ , „ . d ^ 17. ^ 

c \^ Rl[Ep) Rl(Ep) 



(13) 



Vice versa, one can say that photons with time delay be- 
tween ti and t2 are produced by protons of energy Ep only 
if these protons have interacted in the shell ri{Ep) — ro iEp ). 
Here ri{Ep) = r{ti, Ep), and r^ and ti are related by eq. (|l3|). 

We note that the relation above in eq. (^^ allows to esti- 
mate the time duration The of the high-energy GRB signal. 
If 77 <C 1, namely the proton deflection is substantial before 
pN interactions take place, one gets The — Rl{Ep)c~^ ~ 
10^ s (Ep/GeV). In other words, ff one observe photons with 
energy E^ ~ 100 GeV (which are mainly produced by protons 
with Ep ~ 1 TeV) the expected GRB duration at that energy 
is ~ 10 days. If protons do not suffer strong deflection be- 
fore interacting, i.e. r; ^ 1, one obtains The — RL{Ep)c~^r]~^ 
This implies that the high-energy GRB duration may become 
shorter and then more advantageous in a observational per- 
spective. 

As next, let us estimate the differential 7-ray fluence on 
Earth F~,{E-y\ti,tj), during the time interval ti — tj. To this 
aim, we calculate Q~,{E~i;ti,tj) relative to photons with arrival 
time between ti and tj, by substituting ro and R with r^ and 



in eq. (|l^) respectively. So 



Fry^Ej', ti,tj) 



Q.,(E-,(1 + z);ti,tj)At(l + z) 
47tDI{z) 



7 cm ^ GeV 



-1 



where At is the proton injection time and 



c / 1/, 



(14) 



(15) 
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is the GRB luminosity distance, with Ho — 70 Km s" 
Mpc~^ and qo — 0.15. Clearly, the differential 7-ray flux 
on Earth at the time ti can be evaluated by <^~j[E^;ti) = 
F^{EY,U,U + dt)/dt. 

We note that by substituting eq. (|l^) into eq. ( [l4| ) the 
injection time At cancels out, so that it does not affect the 7- 
ray signal on Earth, provided the observational time is > At. 

It is interesting and important in connection with the ob- 
servability at high-energy of GRB events, to study in more 
details the temporal behaviour of the high-energy 7-ray flux 
arriving on Earth. To this purpose, in analogy to the GRB 
soft 7-ray observations, we define the time Tso{E-y) needed to 
receive on Earth the 50% of the total number of photons of 
energy . In Fig. 3 we show Tso as a function of the photon 




2 3 

Log(E,yGeV) 

Fig. 3. The duration time T50 is given as a function of E-^ for 



selected values of the cloud density 



10". lO'-*, 10 



10^ 



energy, for different values of the cloud density n]v . As one can 
see, for each density value, there exist a limiting value of E^ 
above which Tso decreases with increasing the photon energy. 
Moreover, Tso rapidly decreases as njv increases, for a given 
value of E^. From the same figure and the discussion above, 
it is clear that the higher is njv, more easily the source is ob- 
servable at high energy. This is essentially due to two reasons: 
the flux depends approximately linearly on um and Tso is 
shorter as rijv increases. It is also obvious that, in order to ac- 
cumulate a significant amount of 7-ray photons, high-energy 
observations must last substantially longer if the cloud density 
Wat is rather low. 

As we stated in Section 1, it is nowadays widely recognized 
that GRBs are cosmological sources. In this case high-energy 
7-rays suffer attenuation (by electron-positron pair creation) 
due to their interaction with photons of the intergalactic in- 
frared radiation field (IIRF). This effect is taken into account 
through the absorption coefficient t{E^, 2), which can be eval- 
uated through the observational data on the IIRF and adopting 
theoretical models of the galactic spectral energy distribution 
of th e IIRF from stars and d ust re radiation (Stecker & De 
Jager 1997, Salomon & Stecker 1998). In Fig. 4, for illustrative 
purposes, we plot r as a function of E^, for some values of the 




Log(E /GeV) 



Fig. 4. The optical depth r for high-energy 7-ray interactions with 
background low-energy photons is given as a function of the energy 
E~f, for selected values of the source redshift 2 = 0.1, 0.5, 1. 



source redshift z. As one can see, for GRB sources at z > 0.5, 
photons with energy E.^ > 100 GeV are heavily absorbed, while 

for closer GRBs with 2 < 0.1 photons with E^ <1 TeV survive 

absorption. So, the differential 7-ray fluence observed at Earth 
will be 



F®(E-,-ti,tj) = F^{E-,-ti,tj)eyip[-~T(E^,z)] 



(16) 



An analogous relation allows to calculate the differential 7- 
ray flux on Earth {E~,\ t). In Figs. 5a and 5b, <I>® at different 
arrival times is given as a function of the photon energy, for 
GRB sources at 2 = 0.1 and 2 = 0.5, respectively. 

In Fig. 6, for the same values of the source redshift 2, the 
differential fluence on Earth is given in different time intervals. 
As one can see in Figs. 5 and 6, due to the presence of magnetic 
fields, GRBs may last much longer in the high-energy band 
than in the soft 7-ray band. 



5. Ground-based high-energy gamma-ray experiments 

Nowadays, the possibility to detect GRBs in the GeV- TeV en- 
ergy range re lies es sentially on ground-based experiments (see 



e.g. Vernetto 200C). Less constrained in size than space born 



detectors, they can explore high-energy 7-rays detecting the 
secondary particles generated by their interaction in the at- 
mosphere: Cherenkov photons in the case of Cher enkov tele- 
scopes as HEG RA AI ROBICC (PadiUa et al 
(McEnery et al. 1999) and Whipple (Connaug 



1998), Mil agrito 



hton et al. |1997D 



or e in the case of extensive air shower array exper iments as 
the Tibet air shower array (Amenomori et al. 1996 ) and the 
planned detector ARGO (Abbrescia 



1996 



Bacci 



1998) 



The Tibet shower array and the HEGRA AIROBICC ex- 
periment have independently reported significant excesses of 
10-20 TeV 7-rays in coincidence with a few GRBs. 

In particul ar, am ong the 57 GRBs detected by BATSE 
(Meegan et al. 1996) in the field of view of the Tibet array, 



some of them show a significant excess of 10 TeV 7-ray events 
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Fig. 5. Differential high-energy 7-ray fluxes on Earth are given at 
different time U = 1, 10, 100, 1000 s. The GRB source redshift is 
fixed at the value z = 0.1 (a) and z = 0.5 (b). 
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Fig. 6. Differential high-energy 7-ray fluences on Earth are given 
for different time intervals {U,tj) = (0-1), (0-10), (0-100), (0-1000), 
(0,00) s. As in Fig. 5a and 5b, the GRB redshift is fixed to be 2 = 0.1 
(a) and z = 0.5 (b). 



with a timescale of ~ 10 s. The statistical significance was 
estimated to be about 6 a. 

There are four GRBs in the field of view of HEGRA array, 
and the HEGRA group observed 11 7-ray-like events above 
20 TeV in a minute bin coincident with GRB 920925c in 



the GRAN AT/ WATCH GRB catalog (Sazonov et al. |1998| ), 
whereas only 0.93 events are expected as background (the sta- 
tistical significance is therefore 5.4 a). The chance probability, 
taking into account an appropriate trial factor, was estimated 
as 0.3%. Seven out of the 11 7-ray events are clustered within 
22 s, suggesting that the GRB timescale is ~ 10 s. 

More recently, the Cherenkov 



telescope Milagrito (McEnery et al. 1999) detected an excess 
of events within the BATSE error box of GRB 970417a, with 
chance probability 2.8 x 10~^. Since all the 54 GRBs detected 
by BATSE and falling in the Milagrito field of view have been 
examined, the chance probability of observing an excess with 
this significance in any of these bursts is 1.5 x 10"'^. Moreover, 



we note that the energy fluence of GRB 970417a may be ex- 
plained with a source redshift z ~ 0.7 and an is otropi c total 
energy in the TeV range, £^'°y > 10^'^ erg (Totani ^000| ). 

From the above mentioned observations we can conclude 
that GRBs emit also in the GeV-TeV energy range, although 
only a fraction of the GRBs detected by BATSE have been 
observed at high-energy. This is consistent with the fact that 
high-energy 7-rays suffer serious attenuation due to the HRF 
radiation, as discussed in Section 4. 

In the framework of our model discussed in the previous 
Sections, the possibility to detect high-energy GRB emission 
depends on the assumed model parameter values, in particu- 
lar, on the GRB energetics £p, the cloud density njv and the 
magnetic field B. 

A systematic study of the parameter range allowing de- 
tection through several experiments (in operation or planned) 
will be presented elsewhere. However, it is likely that the most 
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suitable instrument to detect GRBs at energies < 100 GeV is 
the GLAST satellite, planne d to be launched by NASA in 2005 
(Gehrels & Michelson Il999|) . 



Here, for illustrative purposes, we refer to the ARGO-YBJ 
(Astrophysical Radiation with Ground-based Observatory at 
YangBaJing) detector, under construction at the YangBaJing 
High Altitude Cosmic Ray Laboratory (Tibet, China). It is 
a cosmic-ray telescope optimized for the detection of small 
size air shower and will be devoted to set several issues in 
cosmic-ray and astroparticle physics including, in particular, 
GRB physics in the e nergy range between ~ 50 GeV and ~10 



TeV (see Bacci et al. 199!; and Vernetto 



1999) 



GRBs are detectable at high-energy if the number 7V^ of air 
showers due to 7-rays from one burst is significant larger than 
the fluctuation a = \/Nb of the background shower number 
A^i, due to cosmic-rays with arrival direction compatible with 
the burst position. So, in order to reduce the background and 
increase the detection sensitivitv, a good angular resolution is 
of major importance (see eqs. ( [l^ ) and (^sj)). 

The angular resolution of the ARGO apparatus ipjo ~ 5" 
(corresponding to e = 0.7) and the effective areas A1^j{E-^) 
and A'^Pf{Ep), for 7-rays and cosmic-rays, respectively, have 
been numerically estimated throu gh Mo irtecarlo simulations 
by the ARGO coUaboration (Surdo 2000[ ) Q. 

The number of events detectable by ARGO in the energy 



range Err 



during an observational time Tabs, is cal- 




culated by the following relations: 



Log(n^, / cm 



Fig. 7. The isotropic GRB energetics C^"" (in units of 10^* erg) 
is given as a function of njvi for different intervals of observational 
time (Tojs = 1, 10, 10^, 10^, 10"* s, from the bottom to the top) 
and GRB redshift 2 = 0.1 (continuous lines), 2 = 0.5 (dashed lines) 
and 2 = 1 (dotted lines). Each curve defines the model parameter 
region above which GRBs should be detectable by ARGO apparatus, 
once the observational time and source redshift are given. 



F^{EY,Q,Tobs)AlfAE^)dE^ 



(17) 



for the GRB signal and 



A/-6 = 27r(l-cosV'7o)robs / ^'^^{Ej,)A'^^f{Ep)dEp (18) 

for the cosmic-ray background. Here, the cosmic-ray flux 
( in units of cm s GeV ^ sr ^) is given by (Honda 
et al. [l995|) 



^Cfl.p , I Ci (Bp + rn^)-'-"^^ e-^ for E^ < 75GeV , . 
* ^^^> - ^ C2 E~^-'" for E, > 75GeV , ^^^> 



where x = -1.871/(0.97 + - ttiI), Ci = 1.085 and 

C2 = 2.103. 

We take an energy range between Emin = 10 GeV and 
Emax ~ 5 TeV, although we have verifled that our results do 
not change substantially if 1 GeV < E,nin < 50 GeV (due to 

the effective area suppression) and Emax > 5 TeV (due to the 
signal and background flux suppression). In Fig. 7, assuming a 
detectabihty threshold condition of 4 standard deviations (i.e. 
Nj > 4y/Nb), we give Sp"" as a function of njv, for different 
intervals of observational time and for GRB redshift z = 0.1 
(continuous lines), z = 0.5 (dashed lines) and z — 1 (dotted 
lines). Each curve defines the model parameter region above 
which a GRB should be detectable by ARGO apparatus, once 
the observational time and redshift are given. The horizontal 



^In the energy range 100 GeV - 5 TeV the effective areas in- 
crease from 2 xlO'^ cm^ to 9 xlO* cm^, for 7-rays, and from 5 
xlO*^ cm^ to 8 xlO* cm^ for CRs. 



line in Fig. 7 marks a conservative (low) value for f"" = 10^* 
erg (corresponding to £p ~ 10^^ An erg). Nevertheless, a GRB 
at 2 = 0.1 produces a detectable high-energy 7-ray signal for 
any cloud density value 10* cm~^ < wat < 10^^ cm~^, while if 
it lies at z — 0.5 a density wat > 5 x 10^ cm~^ is needed. Even 

a GRB at z = 1 could be detectable provided un > 10^^ cm~^ 

and/or S^"" >W^^ erg. 

6. Conclusions 

In this paper we adopted Totani's proposal that GRB sources 
release roughly the same amount of energy £ ~ lO^^AH erg 
and that a substantial fraction of this energy is emitted in the 
form of accelerated protons coming out from the GRB source. 

If there exists a dense enough cloud near or around the 
GRB source, protons, interacting with the nucleons in the 
cloud, give rise to pious and ultimately to a high-energy 7- 
ray signal (in the GeV- TeV energy band) within the error box 
of a GRB detected in the KeV - MeV energy band. 

Due to the presence of magnetic fields (with assumed in- 
tensity _B ~ 1 ^G) around the source, i.e. within the proton 
propagation region before pN interactions are in operation, the 
high-energy 7-ray signal results to be delayed and spread out 
over a longer time with respect to the soft 7-ray signal. 

As it is shown in Fig. 3, the time duration of the 7-ray 
signal in the GeV - TeV energy band crucially depends on the 
cloud density njv and lasts much longer with respect to the 
signal in the soft 7-ray regime. 

By assuming that the accelerated proton injection time At 
is much shorter than the proton interaction time, we have cal- 
culated the 7-ray flux on Earth at different times starting from 
the onset of the high-energy signal and the corresponding flu- 
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ences for difFerent integration times. These quantities depend 
on both the cloud density njv and redshift z of the GRB source. 
In particular, due to the 7-ray absorption from IIFR back- 
ground radiation from stars and dust, the flux (and fluence) of 
the high-energy signal is largely suppressed for source distances 
z>\. 

In Section 5 we have analyzed the possibility to detect a 
GRB in the high-energy band referring, in particular, to the 
performances of the ARGO detector under construction in Ti- 
bet. As it can be seen from Fig. 7, there exists a wide range 
of values for the relevant model parameters and njv which 
might make a GRB detectable in the high-energy band. 

As far as the GRB occurrence rate is concerned, we rely 
on th e analysis of the GRB space density done by Schmidt 



199S who assumed a homogeneous sample of GRB derived from 



BATSE DISCLA data and a wide variety of broken power- 
law GRB luminosity functions. As a result of this analysis, 
one obtains that about 13 GRB/yr are expected to occur at 
redshift z < 0.5 and ~ 120 GRB/yr at 2 < 1. 

However, we note that a GRB with z ~ 0.1 is an extremely 
rare event (less than one GRB/yr). Moreover, it is expected 
that only a fraction of the GRB events meets the conditions to 
be observable on Earth since we do not claim that all bursts 
have £p = lO^'^Afi erg and that accelerated protons interact 
with a very dense medium (with 10*cm~^ < njv < lO^'^cm"^). 

Finally, we note that, although the details of the target 
cloud geometry are not important as far as the high-energy 
7-rays flux estimates, the long-wavelength afterglow crucially 
depends on that geometry. In particular, it is expected that a 
long-wavelength afterglow may accompany the high-energy 7- 
ray production only if the cloud covering factor is low enough. 

Acknowledgements. We would like to thank the ARGO Group 
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